Quantum Dimer Model on the Kagome Lattice: Solvable Dimer Liquid and Ising Gauge 
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We introduce quantum dimer models on lattices made of corner-sharing triangles. These lattices 
includes the kagome lattice and can be defined in arbitrary geometry. They realize fully disordered 
and gapped dimer-liquid phase with topological degeneracy and deconfined fractional excitations, as 
well as solid phases. Using geometrical properties of the lattice, several results are obtained exactly, 
including the full spectrum of a dimer-liquid. These models offer a very natural - and maybe the 
simplest possible - framework to illustrate general concepts such as fractionalization, topological 
order and relation to Z2 gauge theories. 
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, Quantum dimer models 
and Kivelson 



(QDM) were introduced by 
'"^okhsar and Kivelson [Q in the context of resonating 
Rralence-bond (RVB) theories for the high-temperature su- 
4^erconductors [|| . Such models are expected to describe the 
dynamics of singlet bonds (dimer) in quantum disordered 
"pgjm-i antiferromagnets. They can describe two generic 
Qjbases: spin-liquids where the system breaks no symmetry 
^ all and dimer (or valence-bond) crystals where long-range 
■ <stiimer-dimer correlations develops. Recently a genuine liq- 
Jid phase with a finite correlation length was found in a 
JDM by Moessner and Sondhi on the triangular lattice js). 
Ii'ch liquid states have attracted a lot of interest because 
they display both fractional excitations and topological or- 
t^er Q . While fractionalization could play an important role 
some theories of high-temperature superconductors |^ H) , 
C^e topological properties of these liquid states have been 
^iroposed as possible devices to implement quantum bits for 
^Euantum computations 

In this Letter we introduce QDM which realize such a 
f^imer-liquid phase. Due to the geometric properties of lat- 
"^tices made of corner-sharing triangles (the simplest two- 
Ji^nensional example being the kagome lattice [^), these 
]iodels are solved exactly. As for the solvable point of the 
JDM on the square and triangular ^ lattices (see also 
iefs. 1^ and Q), the ground-state is the equal-amplitude 
^aperposition of all dimer coverings in a given topological 
Catectorjps). Such state has been first considered by Suther- 
^nd |llO| (SRK state in the following) and is the prototype of 
: 'Tfesonating valence-bond (RVB) state. However, our model 
raas several important differences with previous analogs : 1) 



ot only the ground state but all excited states wave func- 
tions are known. Elementary excitations are (pairs of) non 
interacting and gapped vortices [TT| (called visons in the re- 
cent literature l^). 2) The model can be solved on any geom- 
etry: torus, discs or spheres. This allows to investigate the 
interplay between topology, ground-state degeneracy and el- 
ementary excitations in an very simple way. 3) Dimer-dimer 
correlations are strictly zero above one lattice spacing. This 
makes this SRK wave function the most possible disordered 
dimer liquid state. 4) This state is inside the liquid phase, it 
is not lying at the phase boundary with a crystalline phase. 




FIG. 1: Medial lattice construction. From a trivalent lattice H 
(full lines) we construct a lattice K whose sites (black dots) are 
centers of the bonds. The sites of K are linked together (dashed 
lines) to form triangular plaquettes. The pseudospin variables 
live on the dual of H (h and h'). 



5) It is known |Q that QDM can be obtained as special lim- 
its of Z2 gauge theories, the gauge variable being the dimer 
number on a bond. Here we show a complete equivalence 
with a Z2 gauge theory. This allows to investigate the con- 
finement transition which goes with a dimer crystallization 
in a simple QDM which exhibits a liquid to solid transition 
accompanied by a vison condensation. 

For all these reasons the model we introduce is more than 
a toy model but the simplest possible RVB liquid. It is a 
"free dimer liquid" point in the short-ranged RVB phase. It 
shows the generic properties which characterize that state of 
matter but almost any quantity can be computed exaclty. It 
is a natural starting point for perturbative expansions toward 
more realistic models. 

Medial lattice construction. — The dimer models described 
in this Letter can be defined on lattices made of corner- 
sharing triangle, constructed in the following way. Let H be 
a trivalent lattice (each site has three neighbors). The lattice 
K, where the dimers live, is the medial lattice of iJ, i.e. the 
sites of K are the midpoints of the bonds of H (Fig. If _ff 
is the honeycomb lattice, K is the kagome lattice (Fig. ^ja). 
In the following, unless mentioned otherwise, we use kagome 
for simplicity, where plaquettes of H are hexagons. We stress 
however that all results can be generalized straightforwardly 
to other lattices (squagome lattice for instance, as well 
as one-dimensional examples). 

Pseudospin representation. — Let us begin with the defini- 
tion of a simple dimer model. For each hexagon /i, we define 
an operator a^ih) as the sum of all possible kinetic energy 



terms involving h only: 

32 

The sum runs over the 32 loops on kagome which enclose a 
single hexagon and around which dimers can be moved (see 
Ref for an explicit list). The shortest loop is the hexagon 
itself, it involves 3 dimers. 4, 5 and 6-dimers moves are also 
possible by including 2, 4 and 6 additional triangles (the loop 
length must be even). The largest loop is the star. For each 
loop a we associate the two ways dimers can be placed along 
that loop: \da{h)) and \da{h)). 

For a given dimer covering \D), all the kinetic opera- 
tors in the sum but one annihilate \D). As a consequence, 
cr^{hY = 1. One can further check that these operators 
flip the pseudospin variables cr^{h) introduced by Elser and 
Zcng |l5| (EZ) to label dimer coverings (Hj. It is im- 
portant to note that cr^ operators depend on the choice of 
a reference \Dq) and are not local, unlike (t^. The (j^{h) 
commute with each other. This is not obvious from Eq. P 
and it is most easily demonstrated in terms of the arrow 
representation that we introduce below. 

Arrow representation. — A correspondence between dimer 
coverings on the kagome lattice and sets or arrows as illus- 
trated in Fig. ||a was introduced by Elser and Zeng . Each 
arrow has two possible directions: it points toward the inte- 
rior of one of the two neighboring triangles. If site i belongs 
to a dimer {i,j) its arrow must point toward the triangle the 
site j belongs to. Consider a triangle without any dimer, this 
arrow rule implies that it will have three outgoing arrows. 
Other triangles will have two incoming arrows and one out- 
going arrow. In other words, the number of outgoing arrows 
is constrained to be odd. 

The number of dimer coverings is 

2^/3+1 where N is the 
number of sites [^3|, K has N sites and N arrows, 

2N/3 triangles and one constraint per triangle. However only 
2A''/3 — 1 constraints are independent because their product 
for all triangles is equal to one. We are left with N/3 + 1 Ising 
degrees of freedom. The existence of this arrow representa- 
tion is a central reason for which QDM considerably simplify 
on these lattices. For example tT^(ft,) translates very simply 
in the arrow representation: it flips the six arrows sitting 
around hexagon h, which clearly conserves the constraint 
for all triangles and commutes from hexagon to hexagon. 

Rokhsar-Kivelson point. — Consider the following Hamil- 
tonian: 

Ho = -TY,a^{h) (2) 

h 

where the sum runs over hexagons (sites of the dual of H 
in general). Although very simple in the pseudospin vari- 
ables, this Hamiltonian is not obviously solvable when writ- 
ten with dimer operators. In the pseudospin variables the 
ground-state is a fully polarized ferromagnet in the x direc- 
tion, which is the sum of all pseudospin configurations in the 
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FIG. 2: a) A dimer covering on the kagome lattice (fat bonds). 
The corresponding arrow representation is displayed. The 
Kagome lattice is the medial lattice of the hexagonal lattice 
(dashed bonds), b) A pair of visons is created by applying to 
the SRK wave-function a factor (—1) for each dimer crossing the 
cut il. It is an exact eigenstate of Eq. ^. 

EZ basis. Back to dimers, this is nothing but the sum 
of all dimer configurations in a given topological sector, that 
is a SRK wave- function [|l|. The ground-state appears to be 
unique in each topological sector, which gives a global 4-fold 
degeneracy on the torus. 

Correlations. — Correlations in the SRK state of the 
kagome lattice are particularly simple: irreducible dimer- 
dimer correlations are strictly zero when their correspond- 
ing triangles do not touch. The arrows on two bonds are 
independent provided they are not involved in a common 
constraint, that is a common triangle. As a result, dimer 
on the kagome lattice are the most possible disordered: they 
are independent above a finite distance. 

Gap. — The whole spectrum is known. The operators 
commute from hexagon to hexagon but physical dimer states 
must satisfy Ylh <^^{h) = 1- This constraint comes from the 
arrow representation since Ylh cr^{h) flips all the arrows twice 
and therefore keeps all dimerizations unchanged. The first 
excited state appears not to be a single but a pair of flipped 
hexagons with energy cost A = 4r. 

Visons. — Despite of the simplicity of the model in the 
pseudospin variables, its excited states are not local when 
expressed with dimer degrees of freedom. A a'^ih) = — 1 
hexagon is a vortex excitation (also called vison Q). Con- 
sider a string which goes from an hexagon a to and hexagon 
b (see Fig. ^) and let Q{a,b) be the operator which mea- 
sures the parity ±1 of the number of dimers crossing that 
string. Q{a, b) commutes with all it^(/i), except for the ends 
of the string: CT^(a)ri(a, 6) — — ri(a, 6)cr^(a). A dimer move 
changes the sign of il(a, b) if and only if the associated loop 
crosses the string an odd number of times, which can only 
be done by surrounding one end of the string. This shows 
that r2(a, b) flips the x component of the pseudospin in a and 
b HI] and Vl{a, b) 1-0) is precisely the excited state of energy 
4r discussed above. Visons appear to be perfectly localized 
in this model. 

Spinons. — One can consider the system with two static 
unpaired sites (spinons or holes). As for others QDM at a 
SRK point I, |, the sum of all dimerizations in a given sec- 
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tor remains an exact eigenstate. The energy turns out to be 
independent of the relative distance between spinous, which 
is a strong indication that spinous would be also deconfined 
if they had kinetic energy. Eventually, we note that taking 
a spinon around a vison will change the wave function by a 
— 1 factor, as expected pi] , 

Visons and topology. — On a closed surface of genus the 
spectrum has a degeneracy given by the number of topo- 
logical sectors 2^^ and excitations are pairs of visons, as al- 
ready mentioned. When the sample has edges, the constraint 
H/i cr^ (h) = 1 is not valid anymore. To handle this case we 
introduce a''{h) which flips the arrows along the edge and 
which restores a^{h) Ylf^ a'^{h) — 1. The excitations are still 
pairs of visons but one vison can be located in the hole js) . 
In this case the gap reduces to 2T. In an cylinder geometry, 
as in Ref. Q for instance, we have two sectors and a dou- 
bly degenerate spectrum. It is interesting to note that this 
dimer liquid has no low-energy edge states. 

Liquid - solid transition. — We consider a new QDM which 
is a generalization of Eq. ^: 

7^1 = Ho -J X! (^"{h)a%h') (3) 

{hji') 

where the second sum runs over pairs of neighboring 
hexagons. The ct^ operators are those defined by EZ, they 
depend on the choice of a reference dimerization \Dq). A 
term a^{h)(7^{h') is local, it is = 1 if the arrow which is in 
between h and h' is in the same position as in the reference 
state and —1 otherwise. Since (j^{h)a^(h') is equivalent to 
il,{h,h'), we see that such a term allows to create, annihi- 
late and move visons in the system. From another point of 
view, the J term in Eq. || counts the number of arrows to 
be flipped to recover the reference state. If J goes to +oo, 
this Hamiltonian obviously selects the reference state as the 
ground-state and dimers are completely frozen. In the pseu- 
dospin language Hi is an Ising ferromagnet in transverse 
field, which displays a second order phase transition at a 
critical value of T separating two phases: a ferromagnetic 
phase with (tr^) > and paramagnetic phase with (cr^) = 0. 
From what we know of the F = and F — > +oo limits we can 
identify the first one with a dimer solid and the second one 
with the liquid. From the Ising point of view the solid phase 
is characterized by (a^) > 0. The two ferromagnetic Ising 
states ((cr^) = ±1) correspond to the same dimer state on 
a closed surface but differ along the edge for open systems. 
The Ising magnetization (cr^(/i)) provides a non-local order 
parameter for the dimer solid, {il{h, oo)), which involves a 
string going to infinity. 

Up to a sign, (7^(/i) is the sum of a creation and anni- 
hilation operators of a vison. The dimer solidification can 
therefore be interpreted as the onset of off-diagonal order 
and macroscopic occupation number (of the zero- momentum 
state) for the visons. This QDM realizes a condensation of 
topological defects ("kinks" or visons |^) at the con- 

finement transition. 



Notice that although we call it a solid, the large- J phase 
does not go with a spatial spontaneous symmetry breaking 
since the J part of Tii depends on an arbitrary reference 
state through the operators; this term acts as an external 
potential which tends to pin the dimers along the reference 
state. 

Z2 gauge theory. — The model of Eq. |^ is the dual of a 
Z2 gauge theory pO| (in its Hamiltonian formulation) where 
the gauge degrees of freedom T^{i) live on the bonds of H 
{i.e. sites of K). By definition T^{i) is the operator which 
flips the arrow at site i and gauge-invariant observable are 
made of products of around closed loops. Spatial gauge 
transformations require the x component, which we define 
with respect to the reference state \Do): T^{i) = 1 if the 
arrow i has the same orientation as in |Z?o) and r^(i) = — 1 
otherwise. For every site of H (every triangle of K) the con- 
straint reads {ii)T'^ {i2)T'^ {i^) = 1 where zi, «2 and 13 are 
the bonds of H emanating from that site. This expresses 
the fact that physical states must be gauge invariant. This 
shows a one-to-one correspondence between physical state of 
the gauge theory and dimer coverings of K and the redun- 
dancy in the gauge theory is solved by the dimer coverings. 
We wish to express Ho with the gauge degrees of freedom. 
Since (h) operator flips all the arrows around h, it becomes 
the plaquette operator: 

a^h) = Y[r%^) (4) 

i 

where the product runs over the bonds of H surrounding h. 
The Ising interactions can be written 

a^{h)a'{h')=T''{i) (5) 

where i is the common link between h and h' (see Fig. 
As a result the QDM of Eq. | translates into the Hamilto- 
nian of a Z2 gauge theory (in continuous time), which is a 
manifestation of the well-known duality between Z2 gauge 
theories and Ising models jl^, Such a Z2 gauge theory 
has a confined and deconfined phases. A classical result is 
that they can be distinguished through the expectation value 
of a gauge-invariant Wilson loop ^ 

W{cu)=<l[T^{^)>=<l[a^^)> (6) 

where the first product runs over a close loop dw which sur- 
rounds the area w. W{lij) changes from a perimeter law 
^ exp(— |9w|) to an area law ^ exp(— |aj|) when going from 
the deconfined at large F to the confined phase at large J. 
The right-hand side of Eq. ^ expresses W{uj) as a correlation 
function for the dimer problem. W{uj) moves dimers along 
the loop duj and the mapping to the gauge theory tells us 
that its long-distance behavior characterizes the liquid and 
frozen phases. 

Liquid - crystal transition. — Hi has the disadvantage of 
depending on a specific reference state. To remedy for that 
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we introduce a QDM which is reference free and restores the 
fuU lattice symmetry: 



^2 = ^0- JX1s^(/i) 



(7) 



where s^{h) = ±1 is diagonal in the dimer basis and counts 
a factor —1 per anticlockwise arrow around the hexagon h. 
In the J — > CX3 limit the system selects (h) = 1 everywhere. 
This can be achieved if the system spontaneously breaks the 
translation symmetry and crystallizes in an ordered pattern 
of six-dimer "stars" described in Ref. |2^. These ordered 
states are degenerate with others (~ 2^^ where L is the linear 
size) obtained by shifting the dimers along any straight line. 

We have looked numerically (diagonalizations up to 54 
kagome sites) at 7i2 and found evidence for a single and 
continuous transition from the liquid to the star crystal at 
J ~ r. This appears as a collapse of the first excitation of 
the liquid, which then transforms into a degenerate ground- 
state of the crystal phase. We claim that the critical point 
is exactly at J = F from a duality argument. From the 
arrow representation it is clear that a^{h) and s^{h') com- 
mute except if h and h' are neighbors, in which case they 
anticommute. Indeed the algebraic relations of and 
are completely symmetric. In particular, we observed nu- 
merically that the ground-state energy is exactly symmetric 
with respect to the exchange of F and J (this is not true 
for all the eigenstates) and the critical point must lie at the 
self-dual point J = T. 

Conclusions. — The kagome lattice has the remarkable 
property that dimer coverings correspond to the physical 
states of an Ising model on the triangular lattice and, by 
duality, of a Z2 gauge theory on the hexagonal one. Ex- 
ploiting this, we introduced QDMs for which exact results 
are derived. In particular, we obtained for the first time the 
full spectrum of a QDM realizing a dimer liquid phase. It 
explicitly realizes fractionalized excitations and topological 
order. Through several models we showed that QDM on 
the kagome lattice are very simple and natural tools to in- 
vestigate the connexions between frustrated magnets, RVB 
physics and spin-charge separation. 
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The QDM discussed in this work do not aim at describing 
the physics of the spin-i kagome antiferromagnet. We are 
currently investigating a different QDM for this later prob- 
lem. 

Dimer coverings are grouped in topological sectors by consid- 
ering their transition graphs, that is the set of loops obtained 
by superposing two coverings on top of each other. Two cov- 
erings with only topologically trivial loops are in the same 
sector. 

EZ pseudospins variables are different from the well-known 
mapping [M. E. Fisher, J. Math. Phys. 7, 1776 (1966)] used 
to solve the Ising model on a planar lattices in terms of 
dimers. 

Up to a global sign (reference dependent) fl{a, b) is equal to 
a"(a)a"(6). 
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